BMC Biology

The rat tapeworm, Hymenolepis diminuta, is an important laboratory model for uncovering molecular processes that underly the success of tapeworms as parasites. Despite its importance, a high-quality reference genome for this species is lacking. Here we present a highly contiguous and effectively complete genome of H. diminuta assembled from PacBio HiFi long-read sequencing data. Our primary assembly consists of 7 scaffolds (N50=29.25 Mb) with total length of 186.53 Mb, has only 7 gaps, and contains 95.7% complete Lophotrochozoan BUSCOs. Our assembly allows us to confirm aspects of Hymenolepis genome organization, such as high repeat content and unusual chromosomal ends, and to show that Hymenolepis genomes encode [~]10,000 genes. Together with annotations of nuclear tRNAs, mtDNA protein coding genes, and mtDNA tRNAs, our assembly currently provides one of the most complete genome resources for a tapeworm species and will enable research on parasitism, animal regeneration, development, and evolution.

Non-obstructive azoospermia (NOA) represents the most severe form of male infertility, severely limiting a patient's prospects for biological fatherhood when surgical retrieval fails. However, the true biological limits of NOA remain obscured by the inherent limitations of conventional gamete recovery protocols: standard centrifugation frequently causes substantial cell loss, masking extremely rare sperm, while surgical interventions are constrained by spatial sampling biases. Here we report SpermSeek, an integrated AI-guided microfluidic platform for real-time, non-destructive isolation of single sperm directly from semen. Operating at scalable throughput (0.36 mL/h), the system achieves 98.3% detection precision and a 95.5% target encapsulation efficiency, suppressing background debris. In a 59-patient NOA cohort, SpermSeek detected morphologically identifiable sperm in 64.4% (38/59) of cases, spanning diverse genetic etiologies, including AZFb/c microdeletions, and severe histopathological phenotypes, such as Sertoli-cell-only syndrome (SCOS). Notably, among a sub-cohort of 41 patients who remained consistently sperm-negative despite prior medical or micro-TESE interventions, our platform identified gametes in 53.7% (22/41) of these cases. Comprehensive safety profiling in healthy human donors and wild-type mice confirmed that processed sperm retain high DNA integrity and epigenomic concordance (r=0.98), supporting transgenerational developmental stability in mice. Furthermore, in a 26-patient validation cohort, SpermSeek recovered rare sperm in 11 cases. Utilizing gametes from a subset (n=5), we demonstrated their capacity to support early human embryogenesis, yielding high-quality cleavage-stage embryos with confirmed genomic euploidy. This work establishes a highly sensitive framework for re-examining the biological limits of human spermatogenesis, laying the foundation to expand autologous reproductive options for patients refractory to conventional retrieval protocols.

Conserved dorsoventral patterning systems have been proposed as evidence for a common evolutionary origin of centralized nervous systems in Bilateria, yet functional evidence outside vertebrates and arthropods remains limited. Here, we investigated the role of pax6 in the annelid Platynereis dumerilii using a mutant carrying a 61 bp deletion in the paired-domain coding region. Loss of pax6 disrupted ventral neuroectodermal patterning at 34 hpf, causing a shift in nk2.2 expression, narrowing of the nk6 domain, and downregulation of pax3/7, while msx expression remained largely unaffected. These early patterning defects were followed by selective neuronal abnormalities at 48 hpf, including displacement of TrpH-positive serotonergic cells and loss of posterior hb9-positive motoneuron domains. By 6 dpf, additional defects were observed in TrpH, ChAT, VAChT, and nk2.2 expression, accompanied by severe disruption of ventral nerve cord morphology and loss of the characteristic rope-ladder architecture. Together, these findings identify pax6 as a key regulator linking dorsoventral progenitor patterning, neuronal subtype specification, and nervous system morphogenesis in Platynereis. Our results provide functional evidence that the conserved dorsoventral patterning network plays an essential role in annelid ventral nerve cord development and support the view that important components of bilaterian nervous system patterning predate the divergence of major animal lineages.

RNA interference (RNAi) is a widely exploited reverse-genetics tool with potential uses for disease control. Successful RNAi has been reported in trematode-vectoring snails, but the composition of RNAi effector-encoding gene complements, a key driver for RNAi efficiency, remain unstudied in these species. Using bioinformatics and comparative genomics, we searched for orthologues of 115 RNAi effector sequences in genomes or transcriptomes of four snail vectors: Biomphalaria glabrata, B. pfeifferi, Bulinus truncatus, and Lymnaea staginalis. Gene expression patterns of selected RNAi effectors were then examined across developmental stages and tissues of the model B. glabrata snail. At least 74 RNAi-related proteins were conserved across all four species, including core components known to be essential for gene silencing. Classical systemic RNAi-deficient (SID) genes that facilitate systemic RNAi in other systems were absent, suggesting that alternative pathways may compensate for dsRNA uptake and transport. Core effectors of secondary RNAi amplification and heritable RNAi were not detected. Expressions of Dicer-1, Argonaute-2, and the exonuclease Eri-1 did not vary significantly with snail size. A putative RNAi-inhibiting Staufen orthologue showed elevated expression in the ovotestis, while another putative cholesterol-interacting gene was overexpressed in the trunk tissue and may partly contribute to RNAi import. Altogether, our results present the most comprehensive overview of RNAi pathway effectors in major intermediate snail hosts for trematodes. The findings underscore the likely broad potential for RNAi use in trematode intermediate hosts as an experimental tool and potential control method.

TULIP (TUbular LIPid binding) domain proteins (TDPs) are found in all living organisms including bacteria. They have various documented functions, some of which clearly related to their intra- or extracellular lipid transfer activities. Extracellular, BPI-related TDPs of insects (B-TDPs, also known as Takeout-related proteins), are often found in chemosensory organs, but little is known regarding their exact location or how they could contribute to olfaction or gustation. We have surveyed and updated the full set of Drosophila B-TDPs and found that roughly 50% are overexpressed in chemosensory organs. Focusing on three genes clustered on the third chromosome, we provide evidence that at least one of the encoded proteins is secreted in the lymph cavity housing the dendrites of olfactory neurons. Biochemical data give support for a putative function of B-TDPs as odorant transporters, but loss-of-function analyses also hint to a potential role as a barrier against plant-emitted terpenoids.

In trypanosmatids genes, transcribed by RNA polymerase II do not have canonical promoters and are organized into directional gene clusters that mature into monocistronic transcripts by a co-transcriptional process known as trans-splicing. Even though gene expression is regulated mainly post-transcriptionally, it is currently understood that chromatin and epigenetics are also involved in this regulation. In eukaryotes, specific signals are normally required for the occurrence of an appropriate transcription initiation. Among them, trimethylation of histone H3 in lysine 4 is the most conserved signal normally detected at transcription start sites of actively transcribed genes. Unlike many model organisms, trypanosomes do not have defined promoters. Instead, transcription initiates in a bidirectional manner from dispersed regions coincident with divergent strand switch regions located between directional gene clusters (DGCs). In T. cruzi, H3K4me3 was observed at the origins of transcription coincident with divergent strand switch regions (dSSRs) in epimastigotes, but it has not been mapped throughout the whole genome at base-pair resolution or in other life stages so far. Here, we set up the CUT&RUN technique for T. cruzi epimastigotes and trypomastigotes. Consistent with a predominant post-transcriptional regulation along the life cycle, we did not find significant differences between life stages. We corroborated that H3K4me3 is enriched at dSSR adjacent to actively expressed DGCs. Moreover, we noticed that this histone mark exhibits different patterns that correlate with the genomic span of the transcription initiation regions and with transcriptional activity. Furthermore, we unveiled that the most actively transcribed DGCs are associated with shorter dSSRs and are located within the core compartment of the genome displaying a more accessible chromatin.

Evolutionary success in mammals requires coordinated regulation of cognitive functions and reproductive capacity. Such coordination must involve shared genes and molecular pathways between the brain and germ cells, yet direct evidence linking cognition to reproduction across species remains limited. Here, proteomic and transcriptomic analyses were performed experimentally in Ovis aries and Rattus norvegicus, while transcriptomic datasets from Mus musculus, Macaca mulatta, and Homo sapiens were analyzed in silico. We identified 8,464 protein-coding genes shared between the brain and sperm/testis and conserved across five species. In rats, 8,444 of these genes were also shared between the brain and the ovary. Functional annotation classified 3,890 genes as associated with both neurological and reproductive functions, and 1,752 as uncharacterized in these contexts, highlighting candidates for future studies on reproductive and neurological disorders. These findings reveal a deeply conserved genetic network linking neurological and reproductive systems, underscoring the evolutionary interplay that supports mammalian fitness.

Microalgal metabolism relies on their chloroplasts, and involves both nucleus and plastidial-encoded proteins of various evolutionary origins. The plastidial ATP synthase complex is a key player in photosynthesis, and has been extensively studied in plants. However, our knowledge in other photosynthetic eukaryotes remains limited, despite their importance in marine environments. Here, we report the characterisation of a novel homologue of the F-type ATP synthase alpha subunit, hereby named xATPA, widespread in microalgae but absent from other photosynthetic organisms. Comparisons of xATPA sequences and predicted structures revealed a specific feature, the bump domain, and highlighted the absence of an ATP-binding site. We assessed xATPA prevalence in microalgae in the global ocean using environmental data from Tara Oceans, with a particular focus on diatoms, and demonstrate that its expression is associated with polar summer conditions. Using a reverse genetic approach in the model diatom Phaeodactylum tricornutum, we show that xATPAP t has a plastidial localisation, and that xATPA KO mutants exhibit growth deficiencies in a combination of low temperature, low salinity and constant light, consistent with environmental analysis. Surprisingly, both RNAseq and physiological assays suggest that xATPA is not involved in ATP synthase functions. On the other hand, xATPA interacts with other F1 ATP synthase subunits in vitro, which we suggest forms transient unassembled complexes. This study hence represents a comprehensive analysis of a novel protein from the environment to the lab, and reveals a new player in the plastidial physiology of eukaryotic microalgae.

The electrophysiological process of auditory transduction in insects remains largely conjecture due to the unknown role of ion channels localised to the cilia, but experimental evidence supports either NompC or Nan-Iav as the auditory mechanotransduction ion channel. Here, we knocked down two key genes that code for the two candidate sound-activated ion channels using dsRNA-mediated RNA interference. We measured sound-evoked activity of the auditory nerve and intracellular electrical currents from the ciliated ending of individual auditory receptors. We found that the sound-evoked nerve activity was reduced in nompC, nan and ift88 knockdown. Using whole-cell patch clamp recordings we found that nompC and nan knockdown resulted in reduced sound-evoked transduction current. Stochastic depolarisations hypothesised to be mediated from one of the candidate mechanotransduction ion channels, either NompC or Nan-Iav, where not affected by knockdown of either channel. The discrete depolarisations are therefore mediated through another unidentified ion channel. We test the hypothesis that discrete depolarisations are graded action potentials that travel toward the soma through noise analysis of the transduction current and analysis of discrete depolarisations to voltage-steps. As a positive control we also knocked down ift88, a protein essential for transporting proteins, including ion channels, along the cilium and found both the transduction current and the discrete depolarisations decreased. Key pointsO_LIInjection of dsRNA decreased RNA of nompC and nan C_LIO_LISound-evoked nerve activity is reduced for RNAi-mediated knockdown of nompC and nan C_LIO_LINompC and Nan both contribute to the transduction current C_LIO_LIThe stochastic discrete depolarisations are not due to NompC or Nan-Iav ion channel but to a third unidentified ion channel. C_LIO_LINoise analysis of the transduction current and the discrete depolarisations suggests they are graded action potentials that travel in the direction of the soma. C_LIO_LIKnockdown of ift88 reduced both the transduction current and discrete depolarisations. C_LI Significance StatementInsects are important to understand, economically, agriculturally and medically. However, we still do not understand fundamental aspects of how insects detect their own body movements, vibrations and sound. These senses are detected by insect chordotonal organs, specialised miniaturised mechanoreceptors that convert movements into electrical signals through specialised ion channels. Previous experimental work has advocated either NompC or Nan-Iav as the mechanosensitive ion channel. Here, for the first time, we reduced the expression of both nompC and nan and measured the sound-evoked transduction current directly from neurons in a specialised auditory chordotonal organ. In contradiction to previous studies, we show that both ion channels contribute to the transduction current and find that electrical signals termed "discrete depolarisations" travel toward the soma.

The haplochromine cichlid fishes of Lake Malawi form one of the most dramatic examples of recent rapid radiation in vertebrates. Here we describe nine new diploid telomere-to-telomere (T2T) genome sequences generated using ultra-long ONT reads, which include 78 ungapped chromosomes. We provide accurate annotations of transposable elements and tandem repeats, identify rDNA cluster regions and putative centromeres, and confirm previously reported large chromosomal inversions. The putative centromeres are primarily composed of satellite tandem arrays of previously reported 237 bp repeats, but notably on most chromosomes these are organised in a novel structure in which four blocks of satellites in alternating orientation are separated by an inverted pair of ~15 kb sequences we term 'centroids', which have similarity to a non-autonomous DNA transposable element and containing potential CENP-B binding boxes. The methylation dip region indicating the likely active centromere always lies between the centroids, whose separation is almost always around 200 kb (interquartile range 151-221kb). A structurally equivalent but non-homologous organisation is seen in the distantly related Etroplus cichlid genera from South Asia. By comparing these structures across chromosomes and species, we suggest how they may have evolved, and potentially how they could contribute to the rampant sympatric speciation seen in these species, based on meiotic drive and chromosome missegregation.

Reporters are widely used in Drosophila genetics to visualize gene expression and cell lineages. However, uncharacterized limitations in specific reporter lines can lead to data misinterpretation. Here, we identify a consistent, driver-independent tdTomato signal in the adult proventriculus from the widely used lexAop-tdTomato.nls reporter line. This signal was observed across multiple lexA driver combinations and was directly detectable in lexAop-tdTomato.nls responder-alone adult proventriculi lacking any lexA driver and without antibody staining. In contrast, no comparable native red fluorescence was detected in larval proventriculi under the same no-antibody imaging condition. Mouse and rabbit anti-RFP immunostaining further supported the presence of proventriculus-associated tdTomato/RFP antigen in adult responder-alone animals. In larval responder-alone proventriculi, antibody-amplified staining was antibody-source-dependent: a detectable signal was observed only with rabbit anti-RFP, whereas mouse and rat anti-RFP produced no reliable detectable signal under the same staining condition. A driver-matched comparison using lexAop-RFP.nls did not reproduce the proventricular signal, arguing against detectable ectopic activity of the tested lexA driver in this tissue. However, because lexAop-tdTomato.nls and lexAop-RFP.nls differ in reporter/transgene architecture and possibly genomic insertion context, the underlying cause cannot be assigned specifically to the lexAop sequence. Our findings highlight the necessity of including driver-negative and no-antibody controls when using this reporter line in adult Drosophila proventriculus and gut studies.

The intestine plays a crucial role in physiology, nutrition, and immune function, but intestinal anatomy and cell types have yet to be fully characterized in many fish species, the most diverse group of vertebrates. To address this gap, we characterized the structure and composition of the intestine of threespine stickleback (Gasterosteus aculeatus), an emerging model teleost in biological research. Using histology, myeloperoxidase staining, single-cell RNA sequencing, and RNA in situ hybridization, we defined major intestinal epithelial, immune, stromal, and stem/progenitor populations. Goblet cells were abundant in proximal and hindgut, while myeloperoxidase-positive granulocytes were evenly distributed throughout the intestine. To facilitate future experimental studies of stickleback intestinal function, we also developed the first intestinal organoid culture from stickleback and show that these cultures recapitulate epithelial architecture and retain expression of canonical intestinal epithelial markers. This organoid platform enables future functional studies of mucosal immunity, host-microbe interactions, and intestinal physiology in stickleback and related teleosts. Together, our integrated approach provides a comprehensive cell atlas and a novel experimental model for studying digestive and immune functions in threespine stickleback.

Marine photosynthetic microbes contribute substantially to global primary production, yet many algal lineages still lack reference genomes with the continuity and annotation quality required for fine-scale structural, regulatory and comparative analyses. Ostreococcus tauri, one of the smallest known free-living photosynthetic eukaryotes, has been a model marine picoeukaryote for over two decades. Despite successive improvements to its historical reference genome, previous assemblies retained hundreds of gaps and incomplete genes, hampering high-resolution genomic analyses. Here, we present O. tauri RCC4221 genome version 2026, a telomere-to-telomere assembly of all 20 chromosomes spanning 13.34 Mb with no gaps. This assembly combines PacBio long-read sequencing, Illumina short-read polishing, correction of unresolved regions guided by independent Nanopore-based assemblies. The updated reference supports a curated annotation comprising 7,683 protein-coding genes, 48 tRNA genes, 3 rRNA operons, 116 medium-sized noncoding RNAs, one signal recognition particle RNA and 138 small nucleolar RNAs. It also improves gene-model integrity and recovers candidate coding loci absent from the 2014 reference. Structural analyses resolved the organization of the two atypical low-GC chromosome 2 and 19 that contain duplicated regions that were collapsed or misrepresented in previous assemblies. Finally, bisulfite sequencing and PacBio SMRT sequencing revealed a dual DNA methylation landscape, with CG-context cytosine methylation concentrated in gene bodies and N6-methyladenosine (m6A) enriched at the start codon. The updated O. tauri 2026 assembly provides a complete and curated reference resource for chromosome biology, comparative genomics, epigenomics and RNA biology in a model marine picoeukaryote.

Cnidae are explosive harpoon-like organelles localized within stinging cells, or cnidocytes, of the phylum Cnidaria (jellyfish, hydroids, sea anemones, and corals). These unique Golgi-derived vesicular structures define the phylum and are prominent examples of an evolutionary cellular novelty. While recent studies have focused on the developmental specification and regulation of cnidocytes more broadly, less is understood about gene expression patterns, structural variations, and toxin repertoires within distinct cnidae subtypes. Here, we determine the transcriptional profile of two major cnidae subtypes in the sea anemone Nematostella vectensis, nematocytes and spirocytes, using the cnidae-specific structural family of proteins called minicollagens. We first define the in vivo expression patterns for three known and three uncharacterized minicollagen orthologs. We show that four minicollagens are broadly expressed throughout ectodermal cnidocytes in developing larvae and primary polyps while two others are restricted to tentacular cnidocytes. Leveraging whole adult scRNA-seq data and two novel transgenic reporter lines, we then demonstrate that the tentacle-restricted cnidocytes are developing spirocytes that are distinguished by expression of the minicollagen NvNcol5. To deepen our analysis of cnidocyte gene expression, we used a customized RNA-FACS-seq pipeline to determine global transcriptional differences between these two subtypes. This approach identified a suite of differentially expressed genes, illuminating spatial and temporal gene expression dynamics across both developing nematocytes and spirocytes. Altogether, our experiments provide fundamental and novel insights into the specialization of cnidarian stinging cells while establishing a rich set of resources for further investigation.

Entamoeba histolytica is a clinically important pathogenic eukaryote and the causative agent of amoebic dysentery. Entamoeba dispar, a nonpathogenic commensal species that resides in the human colon, is the closest sibling species, and serves as an appropriate comparator for genome-wide analysis. Although the genome of E. histolytica is approximately 26.9 Mb, and the largest known genome within the genus, that of E. invadens, is approximately 40.9 Mb, obtaining high-quality assemblies in this genus has remained challenging due to extensive repetitive regions, tRNA gene arrays, and aneuploidy. Here, we used PacBio HiFi sequencing to assemble the genomes of the pathogenic E. histolytica and the nonpathogenic E. dispar. We reconstructed all 36 chromosomes of E. histolytica and 35 chromosomes of E. dispar, assembling each as a single continuous DNA sequence (contig). The two species exhibited high genome-wide nucleotide similarity and conserved synteny at the amino acid level. At one end of each chromosome, we identified tRNA arrays, whereas the opposite end lacked such arrays, resulting in an asymmetric chromosomal architecture. Analysis of unique-read depth revealed widespread aneuploidy in both species: E. histolytica is predominantly tetraploid, whereas E. dispar is diploid, a conclusion further supported by SNP allele-frequency distributions. These assemblies provide a robust foundation for comparative genomics in Entamoeba and offer detailed insights into chromosome-end structure and ploidy.

Chemosensory systems play a central role in host detection, feeding behavior, and habitat selection in hematophagous insects. Here, we performed a comparative evolutionary analysis of chemosensory gene repertoires across 13 species of the Chagas disease vector genus Rhodnius. While gustatory receptors (GRs), ionotropic receptors (IRs), odorant-binding proteins (OBPs), and chemosensory proteins (CSPs) remained globally conserved, odorant receptors (ORs) displayed extensive lineage-specific expansions, tandem duplications, dynamic transcriptomic regulation, and recurrent signatures of positive selection. Major OR expansions were observed in Rhodnius robustus and Rhodnius colombiensis, suggesting increased sensory diversification in ecologically heterogeneous lineages. In contrast, conserved GR1 expression supports the maintenance of ancestral sugar-detection pathways despite hematophagy lifestyle. We further found no evidence of the canonical insect CO2-associated GRs, suggesting alternative molecular mechanisms for CO2 perception in Triatominae. Several receptors, including Orco, also displayed shifts in selective constraints between sylvatic and domiciliary species, consistent with sensory remodeling associated with adaptation to domestic habitats. Together, our results identify ORs as the most evolutionarily dynamic component of the Rhodnius chemosensory repertoire and highlight contrasting evolutionary trajectories among chemosensory gene families during ecological diversification and vector adaptation.

The basidiomycete Ustilago maydis is a well-characterized model organism for studying pathogen-host interactions and of great interest for a broad spectrum of biotechnological applications. We set here to develop light inducible molecular tools to enable dynamic studies on signaling networks and fungi-host communication, and for metabolic engineering approaches. In particular, light-controlled, optogenetic switches provide quantitative, spatio-temporal control capabilities, are minimal invasive and reversible. We engineered two blue light-inducible LOV-domain-based gene expression switches, to up- (Blue-ON) and down-regulate (Blue-OFF) gene expression, and performed a functional characterization in sporidia and hyphae of U. maydis. Profiting from the dynamic control ranges and rapid kinetics, we implemented the optoswitches to control cell morphology by initiating the transition from a haploid sporidial cellular morphotype to filaments upon regulation of the levels of the polarity factor Rac1 and its constitutive active mutant Q61L. In addition to showing how expression level of effectors can be precisely regulated as an approach to understand fungi-plants interaction, we show in two proof-of-principle applications targeted control over U. maydis filamentous fungal invasion of plant tissue and the mechanisms of tumor formation. For this we placed under Blue-ON and Blue-OFF control two U. maydis effectors, See1 (Seedling efficient effector 1) and TIN2 (Tumor inducing 2), and tumor formation was assayed on maize leaves. Taken together, this study established blue-light switches as effective tools to control morphogenesis and pathogenesis in U. maydis.

In most organisms, all cells inherit the same genome, and many mechanisms exist to preserve its integrity across cell divisions. Programmed DNA elimination (PDE), the targeted removal of specific genomic regions from somatic cell lineages during early embryogenesis, is a striking exception. Since its discovery in parasitic nematodes over a century ago, PDE has been observed in diverse eukaryotes, including ciliates, arthropods, and vertebrates. However, the mechanisms, functions, and evolutionary origins of PDE remain poorly understood. Here, we describe the discovery of PDE in three species of the free-living nematode genus Caenorhabditis. Multiple genomic regions are precisely eliminated from somatic cells during early embryogenesis, resulting in chromosome fragmentation and the loss of key germline genes. The sites of elimination are strongly associated with conserved sequence motifs that likely direct DNA breakage. Comparative analyses indicate that PDE was present in the last common ancestor of Caenorhabditis and subsequently lost early during the evolution of many species, including C. elegans. The presence of PDE in the ancestors of one of biology's most important model organisms, together with recent discoveries in other eukaryotic lineages, reveals PDE to be a far more widespread and significant feature of evolution and development than previously recognised.

Black-footed ferrets (Mustela nigripes) are highly susceptible to sylvatic plague caused by Yersinia pestis, but the genetic basis of this vulnerability remains poorly understood. Here, comparative immunogenomic analyses across Carnivora species identified a conserved class of immunoglobulin lambda variable (IGLV) genes with unusually long antigen-binding sites (CDRL1) that are common among Caniformia species but absent in Feliformia species. First discovered in the domestic ferret (Mustela putorius furo), these genes encode tyrosine-rich and anionic motifs resembling the chemokine receptor CCR5 and contain experimentally validated sulfotyrosines previously associated with pathogen-interacting interfaces. Evolutionary analyses revealed distinct selective pressures across Caniformia lineages and showed strong purifying selection acting on long-CDRL1 IGLV genes in mustelids and bears. Antibody repertoire sequencing demonstrated that these genes are actively utilized in expressed repertoires and that their usage correlates with evolutionary conservation. Functional analyses of monoclonal antibodies derived from the long-CDRL1 IGLV gene identified an antibody that significantly reduced intracellular Y. pestis survival in macrophages and revealed a positive correlation between anti-plague activity and sulfotyrosine signal. Notably, all analyzed black-footed ferrets carried a frameshifting deletion in the long-CDRL1 IGLV gene resulting in loss of its expression in antibody repertoires. Together, these findings uncover a germline-encoded immunoglobulin feature conserved across dog-like carnivores and suggest a potential link between antibody germline variation and immune responses to plague.

Meiosis is a fundamental process in producing both sperm and eggs, yet recombination landscapes often exhibit sexual differences, known as heterochiasmy. Since meiotic proteins are generally expressed in both sexes, the molecular mechanism driving heterochiasmy remains elusive. The -kleisin subunit gene of meiotic cohesin, Rec8, is expressed bisexually in mammals, while its putative teleost ortholog, rec8a, is expressed in a female-biased manner, presumably due to the presence of its paralog originating from the teleost-specific whole-genome duplication (TGD). Here, we elucidated the evolutionary history and expression dynamics of -kleisin genes across teleost lineages. Through comprehensive phylogenetic and synteny analyses, we revealed that major teleost lineages retain two copies of rec8 and rad21, with rec8 loci experiencing drastic chromosomal rearrangements immediately after the TGD. Using in situ hybridization and single-cell transcriptome data in medaka and zebrafish, we demonstrated a conserved sexually biased expression pattern: rec8a is predominantly female-biased, whereas rec8b exhibits male-biased expression during gametogenesis. Furthermore, comparative epigenetic analyses revealed that the conserved sexually biased expression is driven by lineage-specific cis-regulatory elements, rather than conserved ones. Motif analyses imply that regulatory rewiring by transcription factors, including foxl2l in particular, might have played a crucial role in the establishment and maintenance of this paralog divergence. Our findings highlight how whole-genome duplication and subsequent genomic and epigenetic rewiring subdivided the bisexual function of rec8, offering insights into sexually distinct meiotic regulation. HighlightsO_LITeleosts possess a unique -kleisin repertoire originating from the TGD. C_LIO_LITeleost rec8 paralogs exhibit conserved sex-biased expression during meiosis. C_LIO_LIDrastic genomic rearrangements after the duplication rewired the teleost rec8 loci. C_LIO_LIThe conserved expression pattern is governed by lineage-specific CREs. C_LIO_LIThose CREs harbor similar types of TFBSs such as Fox-family TFs. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=94 SRC="FIGDIR/small/731870v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@c8c84dorg.highwire.dtl.DTLVardef@1d65668org.highwire.dtl.DTLVardef@c2d732org.highwire.dtl.DTLVardef@1be54a2_HPS_FORMAT_FIGEXP M_FIG C_FIG